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Abstract 



We present results of searches for diphoton resonances produced both inclu- 
sively and also in association with a vector boson {W or Z) using 100 pb _1 of 
pp collisions using the CDF detector. We set upper limits on the product of 
cross section times branching ratio for both pp — > 77+ A 7 " and pp — > jj+W/Z. 
Comparing the inclusive production to the expectations from heavy sgoldsti- 
nos we derive limits on the supersymmetry-breaking scale in the TeV 
range, depending on the sgoldstino mass and the choice of other parameters. 
Also, using a NLO prediction for the associated production of a Higgs bo- 
son with a W or Z boson, we set an upper limit on the branching ratio for 
H — > 77. Finally, we set a lower limit on the mass of a 'bosophilic' Higgs 
boson (e.g. one which couples only to 7, W, and Z bosons with standard 
model couplings) of 82 GeV/c 2 at 95% confidence level. 

PACS number(s): 13.85.Rm, 13.85.Qk, 14.80.-j,14.80.Ly 
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I. INTRODUCTION 



Many processes in extensions of the standard model (SM) result in final-state signatures 
involving two vector gauge bosons, VV + X, where V is either a W, Z, or photon. The 
signature of high mass photon pairs is attractive for searches for new physics as the photon 
is the lightest gauge boson, and hence might be more easily produced in decays of new 
particles. In addition, the photon, being stable, does not decay into many different final 
states as do the W and Z. The dominant SM background process, the production of very 
massive photon pairs (M 77 > 100 GeV/c 2 ), is small compared to the cross-sections for 
producing new strongly-interacting states via quark-antiquark annihilation, making this an 
attractive channel in which to search for new particles or interactions. Examples of possible 
sources of high mass diphoton pairs include sgoldstino production |TJ, interaction terms 
arising from extra spatial dimensions || , a new interaction at a high scale manifesting itself 
as a qq — > 77 contact interaction |3|], a 'bosophilic' Higgs boson or a heavy analog of 

the 7T° that also does not couple to fermions ||. In this paper we focus on the production 
of sgoldstinos and Higgs bosons and their decay into two photons. 

Models with spontaneous breaking of global supersymmetry require a massless and neu- 
tral spin-| particle, the goldstino (G). When gravitation is added and supersymmetry is 
realized locally the gauge particle, the graviton, has a spin-| partner, the gravitino, which 
acquires a mass while the goldstino is absorbed §| . Goldstinos (R-odd) have supersymmetric 
partners called sgoldstinos (R-even) which are expected to be a part of the effective theory at 
the weak scale if gravitinos are very light (< 10~ 3 eV/c 2 ). The simplest model considers two 
neutral spin-0 states: S (CP-even) and P (CP-odd), for which we use the generic symbol </>. 
The mass for these states is completely arbitrary and although initially signals were studied 
in the limit of vanishing masses |K| , we follow the suggestions of Ref . Q and concentrate on 
massive sgoldstinos, = 0(100 GeV/c 2 ). The production of sgoldstinos is dominated by 
the gluon-gluon fusion process while their decay is dominated by two-body decays into a 
pair of gluons, goldstinos, photons, W's, Z's and top quarks. The corresponding branching 
ratios have been calculated |[l l| for two specific choices of parameters, the branching ratio 
into two photons being of the order of a few percent. Limits on the supersymmetry-breaking 



scale v F have been set by the DELPHI Collaboration [|12j for sgoldstino masses up to about 
200 GeV/c 2 . We take advantage here of the higher energy reached at the Tevatron to extend 
the search to much larger masses. 

There are also models in which a Higgs boson could decay into two photons with a 
branching ratio much larger than predicted in the standard model. Figure [I] shows the 
dominant diagrams for production of a standard model Higgs boson (H) in pp collisions. 
The total production cross section is dominated by the gluon-gluon fusion process, and has 
a value of approximately 1 pb for Mn ~ 100 GeV/c 2 |5||n|. Figure ||] shows the dominant 
decay diagrams for a SM Higgs with mass less than ~130 GeV/c 2 . The dominant decay 
mode of the H in this mass range is H — ► bb, with the branching ratio to 77 being on the 
order of 10 -4 . At higher masses, the decays to vector boson pairs WW and ZZ dominate. 
However, some models beyond the standard model introduce anomalous couplings [O] or 



additional Higgs multiplets f|0], enhancing the coupling to photons or suppressing the 
coupling to fermions. The result is a low-mass Higgs boson with significantly increased 
branching ratio to two photons. In the bosophilic models, the coupling to fermions at tree 
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level is set to zero while maintaining the SM coupling to vector bosons. Although the decay 
to two photons proceeds through a higher-order loop diagram, it is the dominant decay for 
M H < M w . For Mu > M w the decay H — > WW* becomes dominant. Since the bosophilic 
Higgs has no coupling to fermions, the gluon-gluon fusion production mechanism is lost and 
the dominant production mode in pp collisions at y/s = 1.8 TeV is associated production 
with a W or Z boson. For Mh = 80 GeV/c 2 , the total associated production cross section 
is about 0.8 pb. The limit set in this paper uses the branching ratios of reference 

Limits on the mass and branching ratios of a bosophilic Higgs boson have been set by 
the OPAL Collaboration assuming SM production of ZH with a lower limit on Mh of 
96.2 GeV/c 2 at 95% confidence level (C.L.) |TJ. More recently, a limit of 100.7 GeV/c 2 at 
95% C.L. |l(| has been set by the ALEPH Collaboration. The DO Collaboration at Fermilab 
has set a lower limit of 78.5 GeV/c 2 at 95% C.L. |J7] in a search at the Tevatron for WH 
and ZH production. 

In this paper we describe a search for departures from SM expectations for both inclusive 
high-mass 77 production and also 77 production in association with slWoyZ boson. This 
search uses 100 ±4 pb -1 of data collected between 1992 and 1995 with the Collider Detector 
at Fermilab (CDF). The photon selection criteria for this analysis were optimized to remain 
efficient for very high energy photons. The analysis is complementary to the previous QCD 
diphoton cross section measurement [T^j. In this present analysis, the photon selection 



criteria have been optimized for high efficiency, taking advantage of the smaller jet fake 
background rate at high Et- The analysis is also complementary to the recent diphoton + 
X search analysis |1{J which searched for non-resonant diphoton signatures, such as ee77$ , T , 
that might arise in gauge-mediated supersymmetric models. 



II. THE CDF DETECTOR 



We briefly describe the CDF detector, which is described in detail elsewhere ||20|| . The 
magnetic spectrometer consists of three tracking devices immersed in the 1.4 T field of a 
3 m-diameter 5 m-long super-conducting solenoid. The magnetic field and three tracking 
devices are all arranged with their principal axis parallel to the proton beam direction (z- 
axis) [plfl. The tracking device closest to the beam line is a four-layer silicon micro-strip 



vertex detector (SVX), used to find secondary vertices, with layers at radii from 2.8 cm 



to 7.9 cm [22]. Surrounding the SVX is a set of time projection chambers (VTX) which 
identifies the pp interaction point (s) along the beam axis with a series of r — z measurements 
out to a radius of 22 cm. The central tracking chamber (CTC) is a 3.5 m-long 84 layer drift 
chamber surrounding the VTX. The CTC wires, ranging in radius from 31.0 cm to 132.5 cm, 
are arranged in 5 superlayers of axial wires alternating with 4 superlayers of stereo wires. 
The calorimeter, which is constructed in projective electromagnetic and hadronic towers, 
consists of the central barrel (\rj\ < 1.1) which surrounds the solenoid, the end-plugs (1.1 < 
\r)\ < 2.4) which form the magnet poles and the forward calorimeters (2.4 < \r)\ < 4.2). Wire 
chambers with cathode strip readout (CES) are located at shower maximum in the central 
electromagnetic calorimeter. These chambers provide a two-dimensional shower profile which 
is used to discriminate on a statistical basis between photons and 7r° backgrounds. Additional 
statistical discrimination is provided by exploiting the difference in conversion probability for 
single photons and pairs from ir° decays in the 1 radiation length of the coil. The presence 
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of a conversion is detected using wire chambers (CPR) located between the coil and the 
central calorimeter. The central muon chambers (\rj\ < 1.1) are located outside the central 
calorimeter to detect particles penetrating the calorimeter. 



III. DIPHOTON EVENT SELECTION 



Photons are identified as a narrow shower in the electromagnetic calorimeter with no 
associated high-P-^ charged particle track. The energy in the hadronic calorimeter and 
adjoining regions of the electromagnetic calorimeter must be small to reject jet backgrounds. 
For high-.&r photons there is a background from 7r° — > 77 decays where both photons are 
very close together. 

The candidate 77 events must pass the diphoton requirements of the three-level CDF 
trigger. The first hardware level requires two central electromagnetic calorimeter trigger 
towers with > 4 GeV. The second hardware level requires two central electromagnetic 
trigger clusters [^] with Et > 16 GeV and a ratio of hadronic to electromagnetic energy 
satisfying Et{HAD)/Et{EM) < 0.125. In the third trigger level, electromagnetic clus- 
ters j24]] are found using the offline reconstruction algorithm and the 16 GeV threshold is 



re-applied to the recalculated transverse energy of the new cluster. 

Offline event selection requires at least two central electromagnetic clusters each satisfy- 
ing the following requirements: Et > 22 GeV, no track pointing at the cluster (or one track 
with Pt < 1 GeV/c), pulse height and cluster shape in the central electromagnetic strip 
chamber (CES) consistent with a photon (to reject 7r°'s and cosmic rays), no additional CES 
cluster in the same 15° azimuthal section of the calorimeter (to reject 7r°'s), and minimal 
energy deposited in the hadronic calorimeter towers behind the cluster. 

Isolation requirements, based on track and calorimeter activity in an rj — cone with 



radius AR = J (A0) 2 + (Ar/) 2 = 0.4 around the cluster, are used to reduce backgrounds 
from jets: EP T (Tracks) < 5.0 GeV/c and (E T (AR < 0.4) - E T (Cluster)) < 2.0 GeV. The 
calorimeter isolation energy is corrected for leakage from the cluster and for pile-up from 
multiple interactions. The efficiency of the calorimeter isolation requirement is studied as a 
function of Et using a sample of electrons from W — > ev events. The efficiency for electrons 
with 30 < E T < 100 GeV is 94.0 ± 0.1% and for electrons with 100 < E T < 200 GeV is 
94.9 ±0.6%. Two requirements reject backgrounds from cosmic rays: there must be at least 
one reconstructed primary vertex within ±60 cm of the center of the interaction region along 
the beam direction, and all energy measured in the central hadronic calorimeters is required 
to be in time with the collision. 

The efficiency to identify a photon passing the above isolation criteria within the fiducial 
region of the central calorimeter is measured using a control sample of electrons from Z° 
decay to be 84 ± 4%. The combined diphoton and event selection efficiency is 63 ± 6% (the 
geometric factor due to the fiducial region is subsumed into the geometric and kinematic 
acceptances, calculated from the Monte Carlo simulation of the detector, as described below). 

Figure ^| shows the invariant mass distribution of the 287 diphoton candidate events that 
pass the selection criteria. A variable bin-width has been chosen to correspond to two times 
the mass resolution (2cr) to enable the observation of narrow structures. 
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IV. BACKGROUNDS 



The dominant backgrounds for this analysis are -y — jet and jet-jet production, where the 
jets have 'faked' photons by fluctuating to a single 7r° or rj, and real photon pairs from 
prompt QCD production. The estimated background from Z° —>■ e + e~ with both electrons 
faking photons is less than 1 event. 

The jet fake rate is measured directly from the data using methods developed for mea- 
surements of the inclusive photon J25[] and di-photon cross-sections [jT8| . For clusters with 
Et < 35 GeV, the lateral shape of the shower in the CES system is used to discriminate 
between prompt photons and photons from 7r° — > 77. Above 35 GeV, where the shapes of 
showers in the CES from photons and 7r°s are indistinguishable, the difference in conversion 
probability of a single photon and a pair of photons (from n° decay) in the material of 
the magnet coil in front of the CPR chambers is used to calculate the single-photon pu- 
rity. These probabilities are used to calculate weights for each event being 'photon-photon', 
'photon-fake' or 'fake-fake'. The result of applying this method to the sample of 287 event 
diphoton candidates is that 183 ±56 ±32 events are 'photon-fake' or 'fake-fake'. This corre- 
sponds to a background fraction of 64 ± 19 ± 11%, where the first uncertainty is statistical 
and the second is systematic. The systematic uncertainty comes primarily from uncertain- 
ties in the modeling of the back-scattering of photons from the electro-magnetic shower in 
the calorimeter into the CPR chambers, and the modeling of the shower shapes in the CES 
chambers. 

The mass spectrum of the jet fakes is determined using a control sample of events enriched 
in fake photons. This sample is made using the same selection requirements as the diphotons 
except that one or both clusters fail the calorimeter isolation requirement. This sample 
contains some real diphotons which fail the isolation requirement. From studies of high- Et 
electrons from W and Z decays, we estimate that 10% of diphoton signal events will end up 
in the non-isolated sample. The mass distribution of the 198 event non-isolated sample is 
normalized to the number of fake events measured in the diphoton candidate sample (183 
events) . 

Two standard model processes make significant contributions to prompt diphoton pro- 
duction: qq — > 77 and gg — > 77. In addition, initial and final state electromagnetic radiation 
from 7— jet production contributes indirectly to the diphoton mass spectrum. In the indirect 
case, several processes contribute to -7 — jet production: qq — > gj, qg — > qj, and qq — > #7. 



These standard model processes are modeled using the Monte Carlo program PYTHIA f26 
with CTEQ4L structure functions |27] and the CDF fast detector simulation. The 77 event 
selection efficiency is determined using the MC and detector simulation, with a correction 
factor of Cmc = 0.76 ± 0.08 applied to account for differences between the detector simula- 
tion and the actual detector performance. These differences are dominated by effects from 
additional low energy tracks from the underlying event and from track reconstruction. The 
correction factor is obtained by comparing the efficiency of the photon selection require- 
ments when applied to electrons from Z° — > e + e~ events from Monte Carlo and data. The 
Z° — > e + e~ events are selected with very loose requirements to minimize any bias in the 
method. A global systematic uncertainty of 13 — 16% applies to these estimates, coming 
from the uncertainty on the correction factor (10%), the modeling of QED radiation (10% 
for diphoton masses below 120 GeV/c 2 and 5% above), the dependence on the structure 
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functions (5%), and the integrated luminosity (4%). 

The total predicted background from fake photons plus QCD diphoton production is 
280 ± 66 events. Figure f| shows a comparison of the diphoton mass spectrum for the 287 
isolated diphoton candidates (points) with background predictions. The shaded distribution 
represents the standard model diphoton prediction from the PYTHIA Monte Carlo program, 
while the unshaded distribution represents the predicted spectrum from jets faking photons. 
The bin-width in this plot corresponds to about 10 times the mass resolution; any narrow- 
width resonance would be seen in the finer binning of Figure |3|. The data are well-modeled 
by the background predictions: above 70 (100) GeV/c 2 we observe 85 (21) events compared 
to a background prediction of 77.1 ± 15.7 (14.7 ± 3.2) events. The numbers of events and 
backgrounds are summarized in Table |. 



V. LIMIT ON INCLUSIVE 77 PRODUCTION 

We first consider the signature of 77 +X. We set limits on the cross section for narrow 
resonances with mass greater than 70 GeV/c 2 [ESJ. The acceptance for diphoton production 



is evaluated using the diphoton decay of a narrow resonance, — > 77, as a model of the 
kinematics for the production and decay of a heavy sgoldstino. The sgoldstino samples are 
generated using the PYTHIA Monte Carlo generator with CTEQ4M structure functions P? 



simulated using the CDF fast detector simulation, and passed through the same event selec- 
tion criteria as the data. The product of efficiency times acceptance increases from 10% at 
75 GeV/c 2 to 16% at 400 GeV/c 2 . The correction factor Cue discussed above is applied 
to the 77 efficiency The acceptance has an additional systematic uncertainty of 4% due to 
the dependence on the structure functions. 

The cross section limit in each mass bin of Table [I] above 70 GeV/c 2 is given by the 
following expression: 

*(PP - 77) < ~ AmJCdt (1) 

where N 95%cl (jj) is the 95% C.L. upper limit on the number of diphoton events in the 
mass bin, e is the selection efficiency, A is the acceptance evaluated in the center of the bin, 
and / Cdt is the integrated luminosity. The upper limit on the number of events in each 
bin is determined using a Monte Carlo technique |29| which convolutes the uncertainties 
(including systematic uncertainties) on acceptance, efficiency and the integrated luminosity 
with the background expectations. The total systematic uncertainty of 12% consists of 
4% from the luminosity measurement, 10% from the selection efficiencies, and 4% from the 
acceptance. Table | provides a summary of the limits. Figure [| shows the cross section limits 
in nine mass bins above 70 GeV/c 2 . For comparison, the cross section times branching ratio 
for pp — > H° + W I Z — > 77 + X production is shown (dashed curve) for bosophilic branching 
ratios ||. The curve corresponding to the standard model branching ratio is not shown, 
being at least one order of magnitude below the bosophilic one. 
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A. Limits on the production of heavy sgoldstinos 



In the scenario in which squarks, sleptons, gluinos, charginos, neutralinos and Higgs 
bosons are sufficiently heavy not to play any role in sgoldstino decays, the most impor- 
tant decays are the two-body decays: <fi — ► GG, 77, gg, jZ, ZZ, W + W~ and //. Three and 
four-body decays are also possible but quite suppressed. Sgoldstino couplings can be pa- 
rameterized in terms of the supersymmetry-breaking scale V~F, the gaugino masses, M3, M2 
and Mi, and a mass parameter, /i a , associated with the charged higgsino. To account for 
the tt(j) coupling for heavy sgoldstinos, two arbitrary free parameters, As and Ap, with the 
dimension of a mass are introduced. We adopt in the following the two sets of choices for 



the parameters adopted in Ref. ||11|| : these choices represent a situation in which sgoldstino 
production is more important than gluino/chargino/neutralino production. The two sets 
correspond to chargino masses of about (220, 380) for case A and about (270, 430) GeV/c 2 
for case B. 

In order for the calculations to be valid, the sgoldstino total width has to be small 
compared to m^. For both parameter sets the decay <fi — > gg dominates, but <fi — > 77 is not 
negligible, being of the order of few percent. 

The dominant mechanism for sgoldstino production is gluon-gluon fusion g + g — > <fi, 
while other associated processes such asq + q^V + (p (V = 7, W, Z) oi q + q — > q + q + <fi 
are suppressed by about four orders of magnitude. The calculation of the production cross 
section has been made at lowest order (LO) however NLO QCD corrections to a(pp — > 
(p) x BR(<p — > 77) are expected to be negligible because they have cancelling effects in the 
cross section and branching ratio. 

Comparing the limits found on the inclusive production cross section to the theoretical 
value of a{pp — > 0) x BR(<p — > 77) bin-by-bin, and considering its 1/F 2 dependence, we 
derive lower limits on yF for sgoldstino masses corresponding to the center of the bin. 
These limits are represented as exclusion regions in the vs y/F space. Figures |6| and [7] 
show these limits for the S- type sgoldstinos. The limits for the P-type (CP-odd) sgoldstino 
are very similar, differing by less than 0.1%. No limit is set in the region > Ms/2, where 
the theoretical calculation may not be valid 



VI. SELECTING 77 + W/Z CANDIDATES 

The inclusive 77 analysis is not sensitive to production of a bosophilic Higgs decaying 
to two photons in the lower- mass region 60-100 GeV/c 2 because the backgrounds from jets 
faking photons and QCD diphoton production are too high (see Figure H). To increase 
sensitivity in this mass region we narrow the signature to be 77 + W/Z. The additional 
requirement of a W or Z boson significantly reduces these backgrounds, allowing access to 
smaller cross sections. 

To achieve a high acceptance for all W and Z decay channels, the vector bosons are 
selected using simple signatures which yield significant background reductions without the 
inefficiency of full reconstruction. The backgrounds from jet fakes and QCD 77 production 
are evaluated using the non-isolated sample of 198 events and PYTHIA Monte Carlo QCD 
background sample used in the inclusive 77 analysis previously described. Backgrounds 
from electroweak processes are found to be insignificant. 
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The vector boson selection consists of the logical OR of three general categories based 
on decay channels as follows: 



1. Central isolated electron (E T > 20 GeV) or muon (P T > 20 GeV/c) for W — > lu and 

z° -> i+r 

2. Two Jets {E T > 15 GeV, \r}\ < 2.0) for W -> gg' and Z° -> gg 

3. $ T > 20 GeV for W Iv and Z° -> i/p 



where $ T is the symbol for missing transverse energy j30 . 

Leptons (e and n) are selected using the isolated central lepton requirements used in the 
'lepton-plus-jets' analysis for the discovery of the top quark fl3~l| . The lepton identification 



efficiencies are measured in data samples of Z bosons decaying to electrons (77.8±0.6%) and 



muons (90.6 ± 0.5%). Jets are identified in the calorimeter using a fixed cone algorithm |32 
with a cone size in r]-(f) space of radius AR = 0.4. Any jet within a radius of 0.4 in r]-<f> space 
of an electron or within a radius of 0.6 of a photon is ignored. Finally the jet-jet invariant 
mass is required to be consistent with a W or Z boson: 40 < Mjj < 130 GeV/c 2 . The 
missing transverse energy is corrected for any high-P^ central muons. Since mismeasured 
jet energies can result in false $ T , events with any jet (E^ > 10 GeV) within 25° of the $ T 
direction are rejected. The same exclusion applies for events with $ T near photons (Ej> > 22 
GeV), electrons {E e T > 20 GeV) and muons (P£ > 20 GeV/c). 

The results of the r y r y+W/Z event selection are summarized in Table [TTT] listing the number 
of events satisfying each W/Z selection requirement. Some properties of the 6 events passing 



the selection requirements are listed in Table [IV] including one event which passes both the 



jet-jet and $ T selection requirements. The highest mass event has a 77 invariant mass of 
137 GeV/c 2 and $ T = 21 GeV. The total estimated background for M 77 > 130 GeV/c 2 is 
0.19 ± 0.12 events, due to standard model 77 production. 



Table |T| also lists the estimated backgrounds from photon fakes, QCD 77 production, 
and electroweak sources, which total 6.4 ± 2.1 events. Fake-photon backgrounds, which are 
estimated from the non-isolated data sample, contribute 1 event to the $ T category and 
3 events to the jet-jet category. Backgrounds from QCD 77, which are estimated using a 
sample generated with the PYTHIA Monte Carlo equivalent to 1 fb _1 of data, contribute 0.8 
events to the $ T category and 1.6 events to the jet-jet category[|. There are small electroweak 
backgrounds, 0.2 ± 0.2 events which contribute to the electron signature from events with a 
W or Z boson produced in association with multiple photons and/or jets. These events only 
contribute in the case where the W(Z) decays to an electron(s) and the charged track(s) 
associated with the electron(s) is not reconstructed. Figure |8] shows the 77 mass distribution 
of events passing all 77 + W/Z selection for the isolated diphoton data and the background 
samples. The mass distribution for the electroweak events is neglected. There is no evidence 
of a 77 resonance in the data. 



There is a small overlap between signatures for the QCD 77 background. 
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VII. LIMITS ON 77 + W/Z PRODUCTION 



We set an upper limit on the cross section times branching ratio for the process pp — > 
77 + W/Z as a function of 77 mass: 

N 95%CL( , w/z) 

°te~rr + w,Z)< t .^ f y» (2) 

where iV 95%CL (77 + W/Z) is the 95% C.L. upper limit on the number of events, e • A is the 
product of efficiency times acceptance, and / Cdt is the integrated luminosity. The number 
of signal events at each mass is taken as the number of isolated diphoton events passing 
the vector (W/Z) selection cuts and falling within a ±3<t(M#) mass window around the 
candidate mass, a being about 2 (3) GeV/c 2 for M H = 100 (150) GeV/c 2 . We calculate 
j\j95%cl a -j. gg^jj m ass, assuming no background subtraction and including a Gaussian sys- 
tematic uncertainty of 15% which includes diphoton selection efficiency (10%), luminosity 
(4%), gluon radiation modeling (11%), and jet energy scale (7%). 

The acceptance is determined from Monte Carlo samples of associated Higgs + W/Z 
generated with PYTHIA and CTEQ4L structure functions [^7j . Figure |9| shows the product 
of the efficiency times acceptance as a function of M# before and after the vector boson se- 
lection cuts. The efficiency times acceptance increases from about 4% for M# = 60 GeV/c 2 
to about 9% for M H > 100 GeV/c 2 . The mass dependence of the acceptance is dominated 
by the photon Et requirement. 

Figure |H] shows the 95% C.L. upper limit on the cross section times branching ratio 
for pp — > 77 + W/Z. The overlayed dashed curve is the prediction for a bosophilic Higgs 
using the branching ratios from reference H and a NLO cross section calculation from 
reference |fL3|| , using the CTEQ4M structure functions ||27|| . A 95% C.L. lower limit on the 



mass of a bosophilic Higgs is set at 82 GeV/c 2 . Table |V| provides a summary of the limit. 

An upper limit on the branching fraction for H — > 77 is obtained by dividing the cross 
section limit on 77 + W/Z by the predicted cross section for W/Z + H production. The 
resulting branching ratio upper limit is shown in Figure O, and lies within the regions 



excluded by OPAL []15|] and ALEPH fllB] . The overlayed dashed and dotted curves are the 
predictions for a bosophilic and Standard Model Higgs boson (scaled up by a factor of 100), 
respectively. 



VIII. CONCLUSIONS 

We have presented results of searches for massive diphoton production both inclusively 
and in association with a high-P^ lepton, $ T , or dijets. The latter channels are sensitive to 
production of a vector boson in association with a Higgs boson which subsequently decays 
to photons. Both the inclusive and exclusive signatures are consistent with predictions from 
standard model sources. In the inclusive channel we set upper limits on the production of 
narrow resonances decaying into two photons. Comparing these limits to a LO calculation 
for massive sgoldstino production we set limits in the range of 1 TeV on the supersymmetry- 
breaking scale \fF for two sets of parameters. In the exclusive channels, we set an upper 
limit on the cross section times branching fraction for pp — » 77 + W/Z between 60 and 
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200 GeV/c 2 . Using a NLO calculation of the SM cross section for pp VH we set a 
95% CL. upper limit on the branching ratio for H — > 77. Between approximately 60 and 
100 GeV/c 2 the upper limit on the branching ratio is less than 1. Using the branching 
ratios of reference || the lower limit on the mass of a bosophilic Higgs is 82 GeV/c 2 at 95% 
CL. 
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FIGURES 




FIG. 1. Diagrams for production of a Higgs Boson in pp collisions: (a) gluon-gluon fusion, 
(b) associated production with a vector boson, (c) and (d) vector boson fusion. In the bosophilic 
models the gluon-gluon fusion diagram is suppressed. 
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(c) (d) 

FIG. 2. Diagrams for the decay of a Higgs Boson to: (a) a quark or lepton pair, (b) vector 
boson pairs (WW/ZZ), (c) via a loop to 77, and (d) via a loop to bb. For a bosophilic Higgs, the 
decay to bb is suppressed relative to 77. 
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FIG. 3. The invariant mass distribution of diphoton candidates (287 events) with a bin-width 
of 4% of the mass. Note that the three highest-mass bins contain one event each. 
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FIG. 4. The diphoton candidate mass distribution is compared with background predictions 
with a bin-width of 20% of the mass. The shaded distribution represents the Monte Carlo prediction 
for QCD diphoton production; the unshaded distribution represents the predicted yield for jets 
faking photons. 
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FIG. 5. Cross section limit at 95% C.L. for high mass diphoton production from a resonant 
state with negligible natural width. The points represent the average mass of the events in each 
bin, but the limits are evaluated at the bin center. The theoretical cross section for a bosophilic 
Higgs boson H is shown as a dashed line. 
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FIG. 6. The exclusion region at the 95% C.L. for S-type (CP-even) sgoldstinos in the Mg-\F 
space for the parameters of set A: M 3 = A s = A P = 400, M 2 = [i a = 300, M x = 200 GeV/c 2 . 
Ms is mass of the S-type sgoldstino. The CDF results are shown as the hatched area; the region 



excluded by results from DELPHI [12] is shown as the solid shaded area. The points represent the 
mass at which the limits are calculated. The boundary T$ = Ms/2 beyond which the model may 
not be valid is also shown. 
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FIG. 7. Exclusion region at the 95% C.L. in the M$-vF space for the parameters of set B: 
M 3 = M 2 = Mi = Ha = A s = A P = 350 GeV/c 2 . The CDF results are shown as the hatched area; 
the region excluded by results from DELPHI p2| is shown as the solid shaded area. The points 
represent the mass at which the limits are calculated. The boundary T$ = Ms/2 beyond which 
the model may not be valid is also shown. 
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FIG. 8. Photon-photon mass distribution compared with background predictions for events 
passing the 77 + W/Z selection. The cross-hatched distribution represents the Monte Carlo pre- 
diction for QCD diphoton production; the shaded one represents the predicted yield from jets 
faking photons. The choppy appearance of the background estimates is the result of low efficiency 
for the W/Z selection. The small electroweak backgrounds are not shown. 
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FIG. 9. Acceptance x efficiency for VH production, with the W and Z bosons decaying via any 
SM decay and the Higgs boson decaying to 77. 
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FIG. 10. Upper limit at 95% CL. on the 77 + W/Z cross section as a function of 77 mass. 
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TABLE I. The number of diphoton events observed, background from jets faking photons, 
'background' from standard model diphoton production, total background, efficiency times accep- 
tance, and 95% C.L. cross section limit for 77 + X production for each mass bin. Mass bins have 
a width of 20% of the bin center. The first two bins are not used for cross section limits, due to 
their low acceptance. 



Set 


M 3 


M 2 


Mi 


Ma 


As 


A P 


units 




A 


400 


300 


200 


300 


400 


400 


GeV/c 2 




B 


350 


350 


350 


350 


350 


350 


GeV/c 2 





TABLE II. The two sets of mass parameters used in the sgoldstino theoretical cross section 
calculations. 
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selection 


Isolated 


Background 




Sample 


Estimate 


Two Isolated Photons > 22 GeV 


287 


280 ± 66 


Central Electron, E T > 20 GeV 


1 


0.2 ±0.2 


Central Muon, P T > 20 GeV/c 








$ T > 20 GeV 


3 


1.8 ±1.3 


2 Jets (E T > 15 GeV, 40 < Mjj < 130 GeV/c 2 ) 


3 


4.6 ±1.9 


Any W/Z signature 


6 


6.4 ±2.1 



TABLE III. Summary of the 77 + W/Z candidate events. The number of 77 candidate events 
passing each of the W/Z selection requirements are listed. There is one event which passes both 
the jet-jet and $ T selection requirements. The background estimates come primarily from fakes 
(non-isolated control sample) plus SM 77 production with a small contribution from electroweak 
sources. Some background events pass more than one of the W/Z selections. 



Run 


Event 


Channel(s) 


M 77 
(GeV/c 2 ) 


Properties 


45219 


277283 


$rjet-jet 


59.1 


$ T = 28.8 GeV, Mjj = 96.1 GeV/c 2 


60597 


119813 




136.8 


$ T = 20.8 GeV 


61514 


9698 


Jet-jet 


48.9 


Mjj = 75.1 GeV/c 2 


68739 


257646 


Electron 


47.1 


E T = 36.1 GeV 


68847 


264160 


Jet-jet 


59.9 


Mjj = 74.6 GeV/c 2 


70019 


155639 


$T 


51.7 


$ T = 22.0 GeV 



TABLE IV. Features of the six events passing the "fj + W/Z selection requirements. The event 
in the electron channel is the ee^^f$ T event J19| |. 
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M 77 


ex A 


(la j dMryry 


(GeV/cr) 




(pb/GeV^J 


60 


0.048±0.002 


1.65 


65 


0.047±0.002 


0.66 


70 


0.055±0.002 


0.57 


75 


0.061±0.002 


0.52 


80 


0.064±0.002 


0.49 


85 


0.066±0.002 


0.47 


90 


0.071±0.002 


0.44 


95 


0.073±0.002 


0.43 


100 


0.074±0.002 


0.42 


120 


0.081±0.002 


0.39 


140 


0.092±0.002 


0.54 


160 


0.087±0.002 


0.36 


180 


0.091±0.002 


0.36 


200 


0.088±0.002 


0.36 



TABLE V. Diphoton mass, efficiency (e) times acceptance (^4), and cross section limit (95% 
C.L.) for associated W/Z + high mass diphoton production. 
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